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We study the behavior of Goldstone bosons in color-flavor-locked (CFL) quark matter at nonzero 
temperature. Chiral symmetry breaking in this phase of cold and dense matter gives rise to pseudo- 
Goldstone bosons, the lightest of these being the charged and neutral kaons K'^ and if". At zero 
temperature, Bose-Einstein condensation of the kaons occurs. Since all fermions are gapped, this 
kaon condensed CFL phase can, for energies below the fermionic energy gap, be described by an 
effective theory for the bosonic modes. We use this effective theory to investigate the melting of the 
condensate: we determine the temperature-dependent kaon masses self-consistently using the two- 
particle irreducible effective action, and we compute the transition temperature for Bose-Einstein 
condensation. Our results are important for studies of transport properties of the kaon condensed 
CFL phase, such as bulk viscosity. 
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I. INTRODUCTION 

Quark matter at sufficiently large densities and low temperatures is in the color-flavor- locked (CFL) state [l| . This 
state is a color superconductor since quarks form Cooper pairs, analogous to electrons in superconducting metals 
[J]. The strong interaction provides an attractive channel in which condensation of quark Cooper pairs takes place. 
Many color-superconducting phases are thinkable, distinguished by their pairing pattern and physical properties. Only 
at asymptotically large densities do we know the ground state of three-flavor quark matter with certainty. The reason 
is that at large densities all quark masses can be neglected with respect to the quark chemical potential /x. In this 
situation, the symmetry group of the Lagrangian is SU{i)c x SU{^)r x SU{3)l x U{1)b with the color gauge group 
SU{3)c, the chiral symmetry group SU{3)fi x SU{3)l and the baryon number conservation group U(1)b- Then, it 
is easy to see that the CFL phase is the only possible color-superconducting phase in which all quasiquarks acquire 
an energy gap in their excitation spectrum. Therefore, the CFL phase has the largest possible condensation energy, 
rendering this phase the ground state at asymptotic densities. 

In the CFL phase, the symmetry group is spontaneously broken to the residual global group SU {3)c+r+l, generated 
by linear combinations of the generators of the original groups. This residual group locks color rotations with flavor 
rotations, hence the name CFL. From the symmetry breaking pattern one concludes that there are 8-1-1 Goldstone 
bosons, an octet originating from the breakdown of the chiral symmetry and a single mode originating from the 
breakdown of U{1)b- The latter is also termed "superfluid mode" since the spontaneous breaking of the particle 
number conservation symmetry accounts for the superfluidity of the CFL phase, just as in the case of more conventional 
supcrfluids like "'He. For inflnitc density and thus a "full" original chiral symmetry, these Goldstone bosons are all 
massless. Therefore, the description of the physical properties of the CFL phase, in particular the transport properties, 
is most conveniently done within an effective theory approach [1, 0, 0] ■ This approach is valid for small energies where 
the excitation of any fermionic degrees of freedom is suppressed due to the energy gap. Hence the CFL phase at zero 
and sufficiently small temperatures can be, for most relevant applications, thought of as a bosonic system of Goldstone 
modes. The octet of Goldstone modes associated with chiral symmetry breaking has the same quantum numbers as 
the corresponding pseudoscalar meson octet in vacuum QCD. Besides the effective fleld theory approach, mesons in 
CFL have also been studied within Nambu- Jona-Lasinio models [1, [^, [13, [HI, [12, [I3l ■ ^^lis approach, which models 
gluon exchange with a point-like interaction, both fermionic quark modes and bosonic meson modes are taken into 
account. 

At moderate densities where the quark masses are no longer negligible, the chiral symmetry is explicitly broken. 
This is certainly the case for densities expected in the interior of neutron stars, where the existence of deconflned 
quark matter is possible [3, [IBl- For these densities, being on the order of several times nuclear ground state density. 
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we expect the quark chemical potential /i to be of order 500MeV. The strange quark mass is expected to be somewhere 
between the current mass ~ lOOMeV and the constituent quark mass ~ 500MeV and can therefore not be neglected. 
Assuming a small explicit symmetry breaking, the Goldstone bosons of the meson octet are still light, however not 
massless. The order of the (zero temperature) meson masses is inverted compared to the QCD vacuum [6|. While the 
superfluid mode remains exactly massless, the next lightest modes are expected to be the neutral and charged kaons, 
K'^ and A'+. In the results of our study, we shall focus exclusively on a system of A'° and K^, one of them being 
condensed. 

The situation in a compact star is further complicated by the constraints of overall electric and color neutrality. 
Together with the heaviness of the strange quark these constraints impose a stress on the CFL phase and eventually 
lead to its breakdown. Whether this breakdown occurs before the transition to hadronic matter is unknown. If so, a 
more exotic color-superconducting phase has to take over. Many possibilities for such a phase have been proposed, 
some of them breaking rotational |16j or translational symmetries. In particular, in the context of Goldstone 
boson condensation, the possibility of a non-zero meson current has been studied (isj . Here we shall not discuss any 
of these possibilities and solely consider CFL quark matter with an isotropic kaon condensate. 

The calculations that we present in this paper are part of an ongoing effort to determine the ground state of quark 
matter at moderate densities using astrophysical observations. To this end, one has to compute the properties of 
candidate color superconductors and test their compatibility with the data. Several approaches have been pursued 
in this direction, focussing on different transport properties of different color superconductors. For instance, cooling 
properties of compact stars have been related to the neutrino emissivity [l9, 20], and glitches (sudden jumps in the 
rotation frequency of the star) have been connected to the rigidity of crystalline color-superconducting phases [2]| . 
Moreover, generic instabilities of so-called r- modes (22j require damping through sufficiently viscous fluids inside the 
stars in order to be consistent with observed rotation frequencies. Therefore, shear viscosity [2^ and bulk viscosity 
[2^ of unpaired and paired quark matter has recently caught attention. Bulk viscosity has been studied for spin-1 
color superconductors [25j as well as for color superconductors where strange quarks remain unpaired (2SC phase) 

For the case of the CFL phase, a calculation has been performed where the kaon modes have been assumed to be 
uncondensed [131 (for a calculation of the bulk viscosity in kaon condensed nuclear matter, see Ref. ^]). This is based 
on the assumption of a sufficiently small kaon chemical potential which, given the uncertainties of this quantity, is a 
possible but rather unlikely scenario. Therefore, it is important to study the opposite case where the kaon chemical 
potential is larger than the kaon mass, allowing for kaon condensation. It is then crucial to have a reliable estimate of 
the critical temperature for this condensation in order to perform calculations for instance of the bulk viscosity that 
are related to a certain temperature regime in the evolution of a compact star. In this paper, we develop a framework 
where both condensed and uncondensed kaons are described at nonzero temperatures and present a calculation of the 
transition temperature. 

The paper is organized as follows. In Sec.|TI]we briefly review the effective theory at zero temperature and discuss its 
simplest extension to nonzero temperature, assuming an ideal gas of mesons. This section mainly serves to introduce 
the notation and to point out the difficulties that arise in the nonzero temperature case. We introduce our self- 
consistent treatment of the temperature-dependent masses in Sec. IIIII for the case of a simple ip'^ theory for a scalar 
field. This serves as the setup for our analysis of mesons in CFL in Sec. IIVI Wc derive the self-consistent equations 
for the and masses and a K'~' condensate. Sec. IIV Al solve them numerically and present an approximate 
analytical expression for the critical temperature. In Sec. |V] we discuss the effect of electric neutrality and the 
presence of electrons on the critical temperature. 

II. BASIC PROPERTIES OF MESONS IN CFL 

We start with explaining how mesons arise in the CFL phase and how they are described in an effective field 
theory. We briefly discuss the zero-temperature predictions of this theory and present the simplest possible nonzero 
temperature approximation. Most of the contents of this section, which introduces the notation used in the subsequent 
sections, can be found in the existing literature. For more details about the following basic properties of mesons in 
CFL, see Refs. 0, 0, 0, H^l, and for more specific studies within the effective field theory see for instance Refs. 
p. 

A. Meson fields from the CFL order parameter 

Diquark condensation takes place in the antisymmetric antitriplet representation [3]c of the color group SU{3)c (we 
neglect the much smaller sextet gaps). Therefore, assuming condensation in the antisymmetric spin singlet channel, 
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we also have to choose the antisymmetric representation in flavor space in order to ensure overall antisymmetry 
of the Cooper pair wave function. More precisely, the order parameter M. lives in the antisymmetric antitriplet 
representation [3]ji © [3]l of the chiral group SU{i)R x SU{3)l, 

Me[3l(E>mR(B[3]L). (1) 

Let us denote the basis for the representations [3]c, [3]^, and [3]l by Ji, and if. respectively (i = 1,2,3), where 
[Jijjk = {Ii')jk = [Ii)jk = —i^ijk- Then, the gauge-variant order parameter for color superconductivity is given by 
a pair Ml = (A^, A^), where A , A^ arc complex 3x3 matrices, 

M = ^fjJ,®lf(B^j:,Jk®Ii- (2) 

In the CFL phase, the order parameter is given by 

A^ = A^ ^ 5,, (3) 

and consequently breaks the group 

G = SU{3)c X SU{3) B, X SU{3) l (4) 

down to the residual group 

H = SU{3)c+R+L ■ (5) 

In other words, the ground state of the system is invariant under transformations of the group H. It is therefore only 
invariant under simultaneous left-handed flavor, right-handed flavor, and color rotations. Hence chiral symmetry is 
broken in the CFL phase. Usually, i.e., in the hadronic phase of QCD, chiral symmetry breaking is induced by a chiral 
condensate of the form {-ipip) . In the CFL phase chiral symmetry breaking is induced by a diquark order parameter 
of the form {i/jiIj). Here the breakdown of the chiral symmetries occurs because of the locking to color degrees of 
freedom. 

Other order parameters than the unit matrix A^ = A^ = 1 are possible. For instance, the so-called 2SC order 
parameter is given by A^ = A^- = SisSj^. In this case, left-handed and right-handed flavor groups are separately 
broken down to SU{2)l and SU(2)j^, respectively. Thus the ground state is still invariant under separate left-handed 
and right-handed flavor rotations. We shall not elaborate on further possible structures of the order parameter. For 
an exhaustive discussion of all possibilities (with the restriction A^ = A-^), sec Ref. [s^ . 

Without explicit symmetry breaking of G, the CFL ground state is degenerate, i.e., it can be rotated by a trans- 
formation of the degeneracy space G/H into an equivalent one with the same physical properties. A (small) explicit 
breaking of G through nonzero quark masses removes this degeneracy, and a rotation with an element oi G/H may 
lead to a physically distinct state, a state where the Goldstone bosons form a Bose-Einstein condensate. Rotations in 
G/H are for example axial color or axial flavor rotations. With V £ SU{3), an axial flavor rotation on = (A^, A^) 
is given by 



F(M) = (A^F,A-^l/-i), (6) 

where A^V^ and A^F^^ are usual matrix multiplications. We see that a state with A-'^ = A^ gets rotated into a 
state with A^ ^ A^. 

One should mention that besides the chiral basis, also a basis of scalar and pseudoscalar order parameters is used 
in the literature [1, H, [l^l- In this basis, the order parameter is a pair (As,Ap), where A''^ = A^ corresponds to a 
vanishing pseudoscalar component, Ap = 0. In other words, for massless quarks one can, due to the degeneracy of the 
ground state, always choose Ap = 0. Now, axial color or flavor transformations generate a nonzero Ap. The explicit 
form of an axial flavor transformation is for instance (As,0) (As cos Q, zA^ sin Q) where V = exp(jQ) G SU{3). 
For the following definition of the meson fields we shall return to the chiral basis and just remark that the definition 
d?]) has the form E = (A^ — Ap)'f(As -I- Ap) in the basis of scalar and pseudoscalar condensates. 

The mesonic field is defined as 

I] = (A-^)tA-^. (7) 

For unitary A^, A^, the 3x3 matrix is unitary, E e U{3). It has 9 degrees of freedom, one of which, the rj', 
corresponds to breaking of the approximate U{1)a symmetry. Since we expect C/(1)a to be explicitly broken at 
moderate densities, we shall ignore the rj' so that E G SU{3). Being a product of a left-handed anti-diquark (or 
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holc-diquark) condensate and a right-handed diquark condensate, S carries quantum numbers of four (anti-)colors 
and four (anti-)flavors. From the above definition of the order parameter we conclude that for example is a 
Cooper pair condensate of two left-handed quarks with colors 2,3 and flavors 1,3 (in an antisymmetric combination), 
is a Cooper pair condensate of two right-handed quarks with colors 2,3 and flavors 1,2 etc. By carrying out the 
matrix multiplication in Eq. ^ we find matrix elements of the form S12 = (AfJ*Af2 + (A2;^)*A^2 + (^3i)*^32 
etc., which can be translated into the respective quantum numbers. Ordering colors as (r, g, b) and flavors as (u, d, s) 
(where we abbreviate red, green, blue by r, g, b, and up, down, strange by u, d, s) we obtain for example 

S32 — > gb^ud]^gb,us]:j + rb,udj^rb,usR + rg,udi^rg,usR^udLUSg. (8) 

We see that the colors cancel each other. This is expected from the transformation properties of E. It transforms 
under the color-flavor group G as S F^EVr where (Vr, Vl) € SU{3)ii x SU{3)l. Hence E is a color singlet. For 
the whole matrix wc flnd the flavor quantum numbers 



/ ds L dsii ds L usu dsL udfj ^ 

usLdsR uslusr usLudR 
\udLdsR udLUSfj udLudu J 



(9) 



We see that the states are composed of four quarks of the structure qqqq as opposed to qq in the case of the usual chiral 
condensate. Nevertheless, we find analogous quantum numbers. For example, the charged kaon states correspond to 
= su, K~ ~ us, generated by the GcU-Mann matrices T4, T5, and the neutral kaon states to K'^ — sd, = ds, 
generated by Tg, T7. 



B. Effective chiral Lagrangian and kaon condensation at zero temperature 

The effective Lagrangian for mesons in CFL is @, 0, [1^, [13] 

C - ^Tr [(9oE + t[A, E])(9oEt ~ t[A, E]t) - i-^a^Ec^.E^] + ^detMTr[A/-i(E + Et)] , 



(10) 



where we have omitted electromagnetic corrections to the charged meson masses [2^ and instanton effects [30|- The 
meson field, defined in Eq. ([7]) can be parametrized as 



E = e^'/f- , 

where 6* is a matrix in the Lie algebra of SU{3). We have defined the "gauge field" 

M2 



A = fiqQ 



2/x 



(11) 



(12) 



Here, Q = diag(2/3, —1/3, —1/3) and M = diag(m„, m^, nis) are the electric charge and quark mass matrices in flavor 
space, and hq the chemical potential associated with electric charge. At asymptotically high density the constants 
Vtt, and a can be determined by matching the effective theory to perturbative QCD [y, |33|, 



J 7T 



21 



! In 2 /x^ 



18 



1 

71 



3A2_ 



(13) 



where A is the fermionic energy gap at zero temperature. Since Lorentz invariance is broken in a medium, spatial 
and temporal components may have different prefactors, ^ I. 

We shall subtract the CFL part of the Lagrangian, which is given by E = 1. Then, we arrive at the new Lagrangian 
(which we also denote by C since from now on this is the only relevant Lagrangian) 



where 



J 7T 



Tr 



Ci = ^Tr 

J 7T 



C — Cq + Ci , 



A^ - {Acose'f - (A sin 6*')^ - 2a{dctM)Ar'^ (1 - cos 61') 



{do cos e'f + {do sin O'f - vl (V cos d'f + (V sin d'f 



« yTr [(^0 cos &) \A, cos &\ 4- (^o sin &) \A, sin B'\ 



(14) 
(15a) 

(15b) 
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with the dimensionless matrix field 9' = 6*//^. Now we separate the zero-momentum mode of the field, 0{K) — > 
0(0) + 9{K) where K = (fco,k) is the four-momentum with the bosonic Matsubara frequencies LOn = ik^ = 2mrT. 
This separation is necessary to take into account Bose condensation as usual {— as in the case of a simple one- 
component ideal Bose gas |3J]). Then, one approximates the zero mode by its expectation value 

9(0) = {9(0)) = $ = cj)aTa , (16) 

where Tq, (a = 1, ... ,8) are the Gell-Mann matrices. In this subsection, we only discuss the zero-temperature case 
for which we only need Cq with 9 replaced by the condensate <f>. In the subsequent sections we shall then turn to the 
remaining terms of Cq that include the fluctuations 9 and the derivative terms in Ci. 

For simplicity, we consider only kaons, i.e., a = 4,5,6,7. Physically, this is motivated by the fact that we expect 
the kaons to be the lightest mesonic degrees of freedom. Moreover, we fix the phase of the condensates by setting 
05 = (t>7 = 0- Then, we denote 0/^+ = (/j^, (j)j^o = (/)g, and 



y^^TT^. (17) 



In this simple case, the matrix $ is has the property <& = $. Consequently, 

(f>\ $<I>v^ 
cos — = 1 - 1 - cos — —5- , sm — = — sm — . (18) 



We may use these identities to compute the zero-temperature free energy density U, which is simply the negative of 
the constant terms in the Lagrangian. After performing the trace in flavor space, we obtain 

/ 2 2 \ '^'k+ , f 2 2 \ 



2 / /■ ^2 ,2 



2 V U 

where we have defined the kaon chemical potentials and masses (squared) 



2 



1-cos^) (^,,+ l^+^^o^) , (19) 



9 2 2 2 

= f^Q + \^ , ^^K<> = , (20a) 

m|.+ = arriiiims + ruu) , m|.o = am„(ms + m^) . (20b) 

We see that these quantities simply arise from evaluating the Lagrangian (jlOp . which, in turn, can be derived from 
symmetry arguments. (Calculations of the meson masses and chemical potentials using an NJL model are in accor- 
dance with these expressions [l3l-) The free energy can be used to find the ground state of the system for arbitrary 
chemical potentials fJ-K+i /^k"- To this end, one has to minimize the free energy through the equations 

^^0. (21) 



and the free energy density is 



By construction, the free energy of the CFL state is given by J7 = 0, corresponding to (f>x+ = (f)xa = 0. If one of the 
condensates vanishes, one of the equations (PT|) is automatically fulfilled. The solution of the other one is 

(h- { 1 for mf > nf 

( for mf>^if 

t/=|_/i(^ fo,„,2<^2 , ^^K\K\ (23) 

By equating the free energies of the two phases 0^+ = 0, 0/^0 ^ and (/)^+ ^ 0, = one finds the condition for 
coexistence of two condensates, 

^^K«(^J'K+ - fn^+f = fJ-K+if^RO - ■ (24) 



This condition can also be obtained by assuming two nonvanishing condensates in Eqs. (j21[) . This results in a phase 
diagram shown in Fig. [1] where we restrict ourselves to fiji+ , /x^o > without loss of generality. Note that we have 
not included any neutrality constraint, i.e., all states in the phase diagram where there is a charged kaon condensate 
are not neutral. We shall use the phase diagram in Fig. [1] as a basis for our nonzero temperature results in Sec. IIVI 
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FIG. 1: Zero-temperature phase diagram for kaon condensation in the jj,i^+-^ifO -plane. No condensation occurs if the chemical 
potential is smaller than the meson mass. Coexistence of the two condensates is only possible along the (solid) line that 
separates the CFL-Jf" from the CFL-A'"'" phase. This line is given by Eq. (|24|l and marks a first order phase transition. For 
large chemical potentials, it approaches the line = fij^o . The (dashed) lines separating either of the two meson condensed 
phases from the pure CFL phase are second order phase transition lines. In the condensed phases, the condensate and free 
energy are given by Eqs. (|22p and (|23|l . respectively. 




C. Expansion in the meson fields and ideal gas approximation 



Before turning to the self-consistent treatment of nonzero temperature effects, let us briefly discuss the expansion 
of the Lagrangian in the fields. Wc shall collect terms up to fourth order in the field. In this subsection, only the 
second order terms shall be evaluated. We shall see that in this simple approximation the system is described by 
decoupled meson components, each corresponding to an ideal Bose gas. Then in Sec. IIVI we include the fourth-order 
terms into our analysis. 

Up to fourth order in the matrix- valued field 6, we obtain from the Lagrangian, given in Eqs. (jl4p and (jlSp . 



Co = ^TT[xe''-iAef] + ^TT 



1, 



:Tr 



{dodf ~ vliWef + 2i{do9)[A,0] 



where we abbreviated 



X = A' 



a{detM)M-^ . 



(25a) 
(25b) 

(26) 



We have neglected all derivative terms of order O'^ such as Ti[{6do9)^] etc. After introducing the condensate by shifting 
the fields as explained above, 9^9 + ^, the Lagrangian becomes 



/: = -t/(<i>) + /:(2)+£(3) + £(4)_ 

Here the tree-level potential for the condensate is 



[/($) = 

The terms quadratic in are 



iTr[X$2 _ (A$)2] + J-Tr 



1(^$)2$2 _ ]_xq>^ - i(A$2)2 
3^ ^ 12 4^ ^ 



(27) 
(28) 
(29) 



with non-derivative and derivative contributions 



(2) 



1 



Tt[X9'^ - {A9j^ 



rTr 



ri, 



[{Aey^^ + {AO, v4$}{6', $} + [A^ye^] 



(2) 



rTr 



{do9Y ~vi{^9Y +2i{doe)[A,9] 



(30a) 
(30b) 
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Here, { — , — } denotes the anticommutator. Wc also find induced cubic interactions, 



£(3) 

and finally quartic terms 



rTr 



^[{Ae, ^ ^Adf{0, $}] - ^X{9^ {0, $}} - ^Ad'A{0, $} 



f 



£(4) ^ ^Tr 
2/^ 



rf 



(31) 



(32) 



-{A0fe-^~-X0^--{A9^r . 

First, we only consider terms quadratic in the fields, i.e., terms oc 9^. This corresponds to the first trace in C/($), 
Eq. jMl), the first trace in C^^\ Eq. ((30a| . and £f \ Eq. ((30b|) . We obtain 



'^('f) = 2i^"''^~^' + ^^'^' + '^2) + «+-M?f+)('/'I + 0i) + «o-/i^o)(0i+</.?) 
H ^ 08 + amuTris \^(p3 - -j= 



amdrris ( 03 + 



(33) 



where fi^+ = IJ.Q + (ml - ml)/{2fi), 



arris (riid + niu). From the terms quadratic in 6 we compute the free 



inverse propagator 5*0 ^, which is an 8 x 8 block diagonal matrix. The charged pion, charged kaon, and neutral kaon 
2x2 blocks read 



[Sq ]ab 



-kl + vlk'^ + ml 



7r+ for a, 6 = 1, 2 
K+ for a, 6 = 4, 5 
A'O for a, 6 = 6, 7 



-2iknfi^ \ 
-kl + i;2fc2 + ml - /if j ' 

The remaining 2x2 block corresponding to 7r°, 77 (a, 6 = 3, 8) is 

-kl + vlk'^ + ams{mu + Wd) 2am.s{md - mu)/^/i 

2ams{md — mu)/V3 —kl + v^k^ + am.s{mu + ?Tic;)/3 + 4am„m(i/3 



(34) 



(35) 



Let us focus on the kaon part of the action. The inverse propagator in Eq. (j34p is obviously the inverse free propagator 
of an ideal Bose gas with chemical potential /i^. Consequently, wc can immediately deduce the thermodynamic 
potential density 17 = —T\nZ, where Z is the partition function. Wc obtain 



where 



e=± • 



(37) 



Here, e = ± accounts for particle and antiparticlc degrees of freedom. In other words, K'^ has a chemical potential 
+IJ,K+ whereas K~ has a chemical potential —^k+ (and analogous for the particle/antiparticle pair K^, K ). Note 
that the zero-temperature part of f2 can also be obtained from expanding the all-order expression in Eq. ()19p . 

Let us comment on the validity of the ideal gas (free particle) approximation. In Fig. [1] we saw that condensation 
at T = takes place for > mf. The ideal gas approximation, however, restricts the chemical potentials to values 
/i? ^ mf . If the chemical potential of any boson becomes larger than its mass then the low-momentum modes have 
unphysical negative energies e^. Therefore, one cannot use the ideal gas approximation for these chemical potentials. 
Note that, although the chemical potential is restricted, the ideal gas can account for arbitrarily large particle number 
densities by increasing the density in the condensate at /x^ = mi. In the physical context of interest, however, the 
mesonic chemical potentials and zero-temperature masses are given and may very well imply /if > mf . This means 
that we cannot use the ideal gas approximation to describe both condensation and nonzero temperature effects. The 
simplest extension of this approximation also fails: suppose we also take into account the second trace in Eq. pOap . 
This leads to corrections to the free propagator 5*0^^. These corrections are proportional to the condensates 0f . For 
small temperatures then the above problem is cured: the resulting bosonic dispersions will be positive for all 
three-momenta k. However, increasing the temperature leads to a decrease of the condensate. Consequently, the 
dispersions approach the ideal gas dispersions for T — > Tc and at sufficiently large temperatures (and still T < Tc) 
negative occupation numbers will occur whenever we consider the case fii > mi. The fundamental problem with these 
approximations is that they are not self-consistent: the excitations of the field around its ground state expectation 
value do not feel the same potential as the expectation value itself feels. To fix this problem wc must set up a more 
elaborate approximation scheme which self-consistently determines the nonzero-temperature meson masses, following 
well-known methods in the context of general bosonic systems [s^, [35| . 
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III. SELF-CONSISTENT BOSON MASSES IN vs* THEORY 

Before we come back to mesons in CFL, we introduce our method with the simple example of f'^ theory for a 
scalar field [13, ■ We shall employ the Cornwall- Jackiw-Tomboulis (C JT) formalism [sl] , as done for example in a 
similar context of a linear sigma model [13] • In contrast to Ref. [13] we have to deal with a fixed chemical potential 
which leads to slightly different self-consistency equations. For a discussion of the linear sigma model with chemical 
potential at zero temperature see for instance Ref. |38l] . Nonvanishing chemical potentials in the C JT formalism have 
been considered at nonzero temperature in Ref. [39|, where a 1/7V expansion in the number of fields TV is employed. 
We discuss the simple case of (p'* theory in detail because we shall reduce the chiral effective theory for mesons in 
CFL essentially to a ip'^ theory, so that we can directly apply the results of this section to the meson case. 



A. Lagrangian and 2PI effective action 



We start from the Lagrangian for a complex scalar field with chemical potential /i, mass to, and coupling constant 

A, 

£ = ](ao - ^^i)p? - \Vp? - m^\^\^ - \\^\^ . (38) 

We introduce real fields via 

ip ^ -^{ipi + iip2) , (39) 

which leads to the Lagrangian 
1 



^ 2 



(40) 

As discussed above for the meson case, we separate the zero-mode (/>i, allowing for Bose condensation, ipi ^ tpi -\- (pi. 
Then the Lagrangian becomes 

£ = -C/(,/.2) + £(2) +£(3)^^(4)^ (4^) 

with 

Um = "^^-^{^l + 4>l) + \{<t>l+ <^lf , (42a) 
£(2) = [-(5o</^i)2-(9o<^2)' + (V</Pi)2 + (V</P2)'-2/i((/.2ao¥'i-¥'i5o(^2) 

+ (to2 - ^_l^){ip\ + ^l) + A(30? + 02)^2 ^ ^(^2 ^ 3^2)^2 ^ 4A0^^2<^^<^2] , (42b) 

£(3) = -A(0i^i+02^2)(¥'? + 'Pi), (42c) 

£(4) = ~\{p\^p?2)\ (42d) 

where we have omitted terms linear in ^pi. These terms vanish upon space-time integration in the action: the field 
(/3i(X) does not contain a Fourier component </?i(if = 0) because this component is separated and given by 
Therefore, the dX-integral over a single field V'i(^) X^x^^o 6xp(i-R' • X)pi(lC) (and an arbitrary number of constant 
factors <^i) vanishes. 

We shall set (p2 to zero which corresponds to choosing a certain direction in the degeneracy space of the ground 
state and thus does not alter the physics. We denote (/> = 0i. The two-particle-irreducible (2PI) effective potential 
V^ff is a functional of the condensate and the full propagator 5, 

Feff [0, S\ = U{cj,) + ^TrlnS-i + h:r[S-\<P)S - I] + V2[<P, S] . (43) 

Here, the trace is taken in the 2-component space of the field {p>i, P2) a-nd in momentum space, and 8^^(41) is the 
tree-level inverse propagator, in this case given in momentum space by the 2x2 matrix 

° I 2ik„fi ^ to2 + A(^2 _ ^2 I ■ ^ > 
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The functional V2[(I),S] is the sum of all 2PI diagrams. We shall use the two-loop approximation for this quantity. 
In this case, V2 [(/ills'] is given by the diagrams shown in Fig. [21 Every line in these diagrams corresponds to a full 
propagator S. The "double bubble diagram" is generated by the quartic term with vertex A, represented by a 
black circle. The diagram on the right, with two black squares, is generated by the terms cubic in the fields ipi, ip2- 
In this case, every vertex contains a condensate cf) as can be seen from £(3) in Eq. (l42cl) . 

The ground state is given by finding the stationary point of Vcs with respect to the variables (j) and S. To this end, 
we define the self-energy 



'ss 



(45) 



In Fig. [21 S is shown diagrammatically. It is obtained by cutting a propagator line in each of the diagrams of V2 [0, S] . 
We shall neglect the right-hand diagram in Fig. [21 which contains two vertices of the order X(f> (black squares). Then 
V2 is a functional only of S, and not of (p. This approximation is particularly simple since the self-energy becomes 
momentum- independent . 



V = 3 





Z = 12 




+ 18 




FIG. 2: Sum of two-particle irreducible diagrams V2[0,S'] in the two-loop approximation and resulting self-energy E. The 
black circle denotes a vertex A while the black squares contain a condensate, \(j>. We shall neglect the second diagram in both 
expressions which originates from the induced cubic interactions. 



B. Stationarity equations 

The self-consistent stationarity equations become 



= 



Sn 



iTr 
2 



(46a) 
(46b) 



The first equation arises from minimization with respect to the condensate, while the second one is the Dyson- 
Schwinger equation, obtained from minimization with respect to the propagator. In order to evaluate these equations, 
we make the following ansatz for the full inverse propagator. 



-K' 



+ Ml - 



-K' 



(47) 



where M_|_ and M_ are masses that have to be determined self-consistcntly. This is the most general ansatz we can 
make, as we prove in Appendix [XI Inversion of yields the propagator 



SiK) 



—2ikQfi 



+ Ml - 



where the excitation energies are given by 



El + T \^^^^El + 5M^ . 



(48) 



(49) 
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We see that the sum and the difference of the two masses (squared) appear in the dispersions, 



, —2 —2 Mi + ^ Mi- 
k^ + M \ M = , = — . (50) 



Let us abbreviate the diagonal elements of the propagator by 
Then, after dividing by (f), the first stationarity equation becomes 



- m2 + + a| ^ [SS+iQ) + S-iQ)] . (52) 
Q 

For the Dyson-Schwinger equation, we need the form of the double-bubble diagram 

K Q 

where the symmetrized vertex tensor follows from the form of the quartic term in the Lagrangian 

£(4) in Eq. (l42dl) . 

Aahcd = J^{5abScd + SacSbd + SadSbc) ■ (54) 

Thus, from its definition (|45p we deduce the self-energy 

S+{Q) + 3S-{Q) ) ■ ^^^> 

Here, the off-diagonal elements vanish due to the antisymmetry of the matrix S. Consequently, with Eqs. and 
(|^7|l the matrix equation (j46bp reduces to the two scalar equations 

Ml = m2 + 3A(^2^A^^[35+(Q) + 5-(g)], (56a) 

Q 

Mi = m^ + X^'' + X^J2[S+iQ) + 3S~{Q)]. (56b) 

Q 

The momentum sum from Eq. (|56aP appears also in the first stationarity equation (|52p . We can use this fact to 
simplify the latter. Moreover, we may write the equations in terms of the mass sum and difference appearing in the 
excitation energies (|49p . To this end, we add/subtract Eqs. (|56ap and (|56bp to/from each other. Finally, we perform 
the Matsubara sum and take the thermodynamic limit to arrive at the set of equations (for details see Appendix |B]) 

+ 5M^ - = , (57a) 

17^ = + 2X(j)^ + 2X I , (57b) 

with the momentum integrals 

r^Tfjpj^Mljl-, , V ). ,;.V./jg^Zlf|) , ' . (58, 



where f(x) = l/[exp{x/T) — 1] is the Bose distribution function. Expanding in the coupling A, we may neglect the 
term A J in the denominator on the right-hand side of Eq. (|57cp . We also approximate 

'-^J (27r)3 E, ■ ('') 
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We have confirmed numerically that, for the results we present below, the exact result is practically indistinguishable 
from the approximate one. We shall see that in the case of kaons in CFL, the quantity that plays the role of the 
coupling A is sufficiently small to use this approximation too. This approximation also has the advantage that the 
Goldstone mode is explicitly gapless. This can be seen by inserting Eq. (|57ap and SA-P ~ XcjP from Eq. (|57cp into 
the excitation energy (|49l) . One obtains e^^g = for all temperatures below Tc- Strictly speaking, Eqs. fST]) violate 
the Goldstone theorem by producing a small energy gap for the Goldstone mode. This shortcoming of the present 
formalism has been observed in previous works in the context of linear sigma models (stI. |40|. (For a self-consistent, 
gapless solution at nonzero temperatures see Ref. [HI.) 



C. Evaluation of stationarity equations 



At zero temperature, / = J = 0, and Eqs. ([57]) yield 

02 = ^^^-^, M' = 2/i2-m2, 5AP^i?-m\ (60) 
A 

Nonzero temperatures have to be studied numerically. In Fig. [3] we present the results of a numerical evaluation 
for the parameters A = 0.01 and ^ ~ l.bm. The mass M decreases with increasing temperature until it equals the 

chemical potential at the critical point. This is expected from Eqs. (|57p . which yield 5M = and M = at T = T^- 
This is also expected by comparison with the theory without chemical potential. In that case we know that the field 
is massless at the critical point. Here, the zero point is replaced by the chemical potential. At the critical point, 
M is non-analytic, indicating the second order phase transition. Above Tc it increases and becomes linear in T for 
sufficiently large temperatures. For T > Eq. (|57aP is automatically fulfilled (we have divided by to obtain this 
equation). Hence, above Tc we simply have to solve Eq. (|57b[) for the single variable M. It is worth pointing out 
that M > fi for all temperatures. Thus, the dispersions are positive for all momenta. Remember that the simplest 
approximations fail to fulfill this elementary requirement, see discussion below Eq. ([37]). The mass SM vanishes above 
Tc since it is proportional to the condensate (f>. 

In the high-temperature limit we may obtain an analytic expression for T^. To this end, we approximate the integral 
/ in Eq. (|57bp by / ~ T^/6. Then, the critical temperature becomes 

Tc' = 3 ^ . (61) 

The minimum of the two excitation branches e^^o shown in the right panel of Fig. [3) We see that one of the 
branches is (approximately) massless, as expected. The other branch has an energy gap of at least 2/i. In a theory 
without chemical potential we would expect the second branch to be massless too at the critical temperature. Here, 
this corresponds to an energy gap of 2/z at T = Tc. 



IV. SELF-CONSISTENT KAON MASSES IN CFL MATTER AND CRITICAL TEMPERATURE 



We are now prepared to turn to the nonzero-temperature analysis of a system of neutral and charged kaons in CFL 
quark matter, taking into account meson condensation. We shall restrict ourselves to the case of only one condensed 
component. This is the most likely scenario in the case of fixed meson chemical potentials. Only after including 
neutrality conditions and neutrino trapping might we expect a coexistence of two mesonic condensates [2^ . We shall 
not consider this possibility in this paper. 

We shall apply the self-consistent method from Sec. IIIII From Sec. |ll] we extract the relevant terms from the 
effective chiral Lagrangian that reduce this theory essentially to a ip* theory. We shall see that basic parameters of 
the chiral effective theory such as the meson masses and chemical potentials play, in certain combinations, the role 
of the coupling constant A of (/s^ theory. The most important difference to (^^ theory is the presence of a second, 
uncondensed boson. Technically, this does not create any additional difficulties. Therefore, we can directly set up the 
self-consistent set of equations upon referring to the previous subsection. Then, we shall solve them and discuss the 
physical results, in particular the value of the critical temperature. 
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T/Tc T/Tc 



FIG. 3: Results of the self-consistency equations (|57|l for theory for a coupling constant A = 0.01 and a chemical potential 
H = 1.5m. Left panel; self-consistent mass M and SM ~ A(/!>^ (solid lines) in units of the chemical potential as functions of 
temperature. The dashed line marks, for comparison, the chemical potential. Right panel; energy gaps e^^g in units of the 
chemical potential. The Goldstone mode is gapless in the symmetry-broken phase, while the other mode assumes its minimum 
2/i at the critical point. 



A. Derivation of self-consistency equations 

We expand the potential U from Eq. (fT9|) up to fourth order in the condensates, 

Uih,^.) + ^^^l + ^^t + ^4 + ■ (62) 

For notational convenience, here and in the following we use the subscript 1 for and 2 for We have abbreviated 



n — ^ lA 2 2\ _ 1 / 2 , 2 2 2 , r \ f l-l-l - fJ-2\^ ,„„^ 



These dimensionless quantities play the role of coupling constants. They can be positive or negative, depending on the 
bare meson masses and chemical potentials . If the condensates were allowed to become arbitrarily large, boundedness 
of the potential (p^ from below would require /3i,/32 > and a > — (/3i/32)^^^. However, here we consider only small 
condensates, since we have employed an expansion of the full potential (the full potential is bounded for all 
possible parameters). Therefore, we do not have any restrictions for the coupling constants. 

In order to determine the 4x4 tree-level propagator for charged and neutral kaons, we consider terms quadratic 
in the field 0, Eqs. ([50]) . We perform the trace in these terms, set the fsr+ condensate to zero, and denote the 
condensate by 4>. Then, the inverse tree-level propagator is block diagonal. 



^02 



where 



^01 - I " „. , >„2^„2j.2^_2_ „2^„ .2 ' (65a) 



°^ ^ [ 2z/i2fco ^k^o+^k^ + ^l^f^l + M' J ' 

For = we recover the ideal gas approximation of the propagator, Eq. (|34p . Upon comparing the inverse propagator 
of the condensed branch, Sq2, with the one from ip^ theory, Eq. (|^ . we see that /?2 plays the role of the coupling 
constant A. Note that the condensate also gives a correction to the bare mass of the uncondensed branch. In this 
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branch, a assumes the role of A (with the important difference, however, that the correction is the same for both 
degrees of freedom) . 

Analogous to above, we use the following ansatz for the full inverse propagator, 

5-'.(^j\°,), (66) 

where 

We have introduced the four self-consistent masses Mi^±, M2,±. Not surprisingly, we shall see below that for the 
uncondensed branch both masses in fact coincide, Mi^^ = Mi^_. Inversion yields the propagator 

51 

'2 



0^ ^J, (68) 



where 

We have defined the excitation energies 



(69) 



ef EE ^Ef + T ^J^^^^Ef + 6Mf , (70) 

with 

/ _2 Mh+Mf ^ Mf.-Mf 

E,^^vlk^ + Ml '■+^ 5Mf^ (71) 

For 5M = 0, one recovers the form of the excitation energy ef — Ei ^ fii. In the following, we shall abbreviate the 
diagonal elements of the kaon propagator by 



C± -V ■ -71- ■ - - r -i (70\ 



-k^ + vie + Mi^-f, j 



For the stationarity equation ()46ap we set 0i = and 02 = </> in the potential U from Eq. (I62p to obtain for nonzero (f> 

fij = ml + /320' + T^Y. H^t + + ^2(3^2+ + ^2^)] • (73) 
Q 

Next, we need the self-energy E. The double-bubble diagram has the same form as in Eq. (|53p (now of course with 
a, 6, c, d running from 4 through 7) , where the symmetrized vertex tensor is extracted from the quartic term in the 
Lagrangian, Eq. ((32|) . 



^iabcd — 4! 9/2 ' 

The self-energy is 



^ATaATbT^Td - ^XTaTbTcTd - -^ATaT^AT^Td 



all permutations of a, 6, c, c? > . (74) 



E„fc = 12Aafccd^5cd(g). (75) 
Q 



Since the propagator S is antisymmetric, the off-diagonal elements do not contribute in the contraction with the 
totally symmetric t' 
contractions yields 



totally symmetric tensor h.abcd- Consequently, E only contains the diagonal elements Sf. Performing the traces and 



E ^ diag(E+ S-, S+, S-) = I ^(A+5+ + \-S- + + X'S') , (76) 
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with 

a;,*" = diag(3/?i,/3i,a, a) , A;," = diag(/3i, 3/3i, a, a) , (77a) 
A^ = diag(a, a, 3/32,/32) , Ag" = diag(Q;, a, /?2, 3/32) ■ ('''7b) 

These diagonal matrices yield the effective coupling constants for the respective contributions to the self-energy. We 
illustrate these contributions in Fig. |4l An interaction between kaons of the same species is given by the effective 
coupling /3i, while an interaction between kaons of different species is given by a. Condensation of species i implies 
/ii > rrii and thus a positive "self-coupling" Pi. The self-coupling of the uncondcnscd species, however, as well as the 
coupling a can become negative for certain values of the bare masses and chemical potentials. 




FIG. 4: (Color online) Kaon self-energy S^^j (meaning or S^) according to Eq. (|76[) . The different contributions originate 
from interactions with the same kaon species (red propagators, same subscript) and with the opposite species (blue propagators, 
opposite subscript). The effective couplings for these interactions are /3i (interaction between charged kaons), P2 (interaction 
between neutral kaons) and a (interaction between charged and neutral kaon), see definition (|63p . Moreover, in the case of 
self-coupling, there is a difference between coupling to the same degree of freedom (solid propagators) and the orthogonal degree 
of freedom (dashed propagators), the former being larger by a factor 3 than the latter. 



With the self-energy ([76)) . the Dyson-Schwinger equation (j46bp becomes the set of four equations 

Q 

Ml_ = ml + + ^ ^ [a{S+ + S^) + pi{S+ + S^f)] , (78b) 

Q 

Ml+ = ml + 3/320' + I 5] [a{St + S^) + ^{^8+ -f S^)] , (78c) 

Q 

Ml_ = ml+M^ + ^Y. H^i P2 {S+ + )] . (78d) 

Q 

We now insert Eq. (|78cp into the first stationarity condition ([75)1 . Furthermore, we add/subtract Eqs. (|78cp and (|78d[) 
to/from each other. Finally, from subtracting Eq. (|78b[) from Eq. (|78ap we conclude Mi_+ = Mi__ = Mi such that 
we can drop Eq. (j78bp . Then, after performing the Matsubara sum, we arrive at the following set of four equations 
for the four variables 0, Afi, M2, SM2, 

III - = /32(/)2 , (79a) 

mI = ml + aif + al2 + 2l3ili, (79b) 

mI = m' + 2/32</'2 + a/i + 2/32/2 , (79c) 
5Ml = /320^ (79d) 

with the momentum integrals 

e=±' 



/.4EKW) + sr«j)i.i:/|^/M)^ (80) 



We have made use of the same approximations as explained below Eq. (|58p in the context of Lp^ theory. In this 
approximation, one of the branches of the condensed species is explicitly gaplcss, e;^(fc = 0) = for all temperatures 
smaller than the critical temperature. 
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B. Solutions for zero and nonzero temperature: melting the kaon condensate 



We may now solve the system of equations (j79p . In these equations, the kaon masses and chemical potentials are 
kept as parameters. A numerical evaluation should be done for values of these parameters that are realistic for the 
interior of a compact star. However, both masses and chemical potentials are known only with significant uncertainties. 
These uncertainties originate partly from the matching calculations that determine the decay constant /,r and the 
constant a, see Eq. (jl3p . These calculations have been performed using perturbative methods within QCD. Hence 
they are reliable only for quark chemical potentials far beyond realistic values. Extrapolations from these high-density 
values yield j-n ~ 100 MeV and a ~ 0.01 for [i ~ 500 MeV and A ~ 30 MeV. The next source of uncertainties are the 
quark masses that enter the expressions of the meson masses and chemical potentials, see Eq. (j20p . Using a strange 
mass of ms — 150 MeV, light quark masses rrid — 7 MeV, niu — 4 MeV, and /ig = we find that the kaon chemical 
potentials are of the order of tens of MeV, ~ ^2 — 20 MeV, while their zero-temperature masses are of the order 
of a few MeV, toi ~ 5 MeV, TO2 — 4 MeV. Due to the slightly larger charged kaon mass, it is reasonable to assume 
condensation of the neutral kaon, see zero-tcmpcraturc phase diagram in Fig. [TJ We have also neglected instanton 
effects which yield an additional correction to the masses. Its value, however, is poorly known, since high-density 
calculations yield results that strongly depend on the QCD scale Aqcd- The value of these corrections has been 
estimated to be of the order of 10 MeV [30|. 

In our numerical results, we shall set the values of the kaon masses to the above values. As for the chemical 
potentials, we have to recall our expansion in the condensate. This expansion is valid for small condensates and thus 
for small values of 1 — ml/^l. In other words, the chemical potential has to be larger than the mass (otherwise 
there is no condensation) - this is likely to be the case judging from the above estimates. However, it should be only 
slightly larger (otherwise our approximation is not reliable) - this is questionable but not excluded from the above 
estimates. Further studies which take into account all orders in the condensate are required to confirm our estimates 
for the cases where 1 — rnl/fil is not small. Finally, we have to recall that the effective theory itself is only applicable 
for energies smaller than the fermionic energy gap A. Also, the critical temperature T^^^ for the transition from 
CFL to unpaired quark matter of course sets an upper limit for our treatment. This temperature is of the order of 
the zero-temperature gap, T^^^ = 0.57 • 2^/'^A ~ 0.7A according to BCS theory (gauge field fiuctuations, however, 
increase the critical temperature (4^.) Hence, when we find critical temperatures Tc for kaon condensation on the 
order of or larger than A and T^^^, the conclusion is "the critical temperature is in a regime where we do not trust 
the effective theory" rather than "the critical temperature is Tc". These are the premises under which our results 
have to be understood. 

Let us first discuss the simple case of zero temperature, where li in Eqs. (|79p vanishes. In this case, 

(81a) 
(81b) 











2,^(2 2\ 
P2 




mI 


= 2/^2 - m\ , 






2 2 

= M2 - ™2 • 





(81c) 
(81d) 

In Eq. ((8Ta|) we have expanded the condensate (\? with respect to the smallness parameter 1 — 7712/^2 ^^d used the 
definitions for a and /?2 from Eq. ([63|) . The result is in agreement with the expansion of the exact result ((22|) . We 
may compare the zero-temperature masses with the bare masses. For the condensed meson the self-consistent mass is 
larger, M2 > m2. For the uncondensed meson, this is not necessarily true since the correction through the condensate 
in Eq. (j81bp becomes negative if a < 0. We also see that the uncondensed kaon becomes gapless in the case of equal 
masses and chemical potentials, ?7ii = 7712, /ii = ^2- In this case, we conclude from Eq. (j63p that a = /3i = /?2 and 

thus m1 = fif. Consequently, e^(A; = 0) = 0. This reflects the fact that exact isospin symmetry, rriu — md and thus 
7771 = m-2 and fii = ^2, together with (neutral) kaon condensation leads to two exact Goldstone modes [43| . 

Let us now turn to nonzero temperatures. Before solving the equations for arbitrary temperatures, we may first 
derive an approximate expression for the critical temperature. We do so for temperatures large compared to the kaon 
masses and chemical potentials. Wc set = in Eqs. ([75]) and approximate 



Then, we obtain for the critical temperature 



a + 2/32 5(/ii - ml) + {fij - mj) + 6^2(^1 + 2^2) 
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where we used the definitions of a and /32 in Eq. (|63p . The masses at T = Tc are 

2/3i 



mI 



2/32 



m: 



(84) 



As it turns out, these high-temperature expressions are a very good approximation to our numerical results. For the 
critical temperature results we shall present below, the approximation is practically indistinguishable from the full 
numerical evaluation. From Eq. (|83p we see that for a + 2p2 < there is no solution for the critical temperature (of 
course assuming that Tc is real). We comment on this apparently unphysical result in Appendix ICl 
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FIG. 5: (Color online) Self-consistent kaon masses and energy gaps as a function of temperature for mi — 5 MeV, m2 — 4 
MeV, fj,i — j-L2 = 4.5 MeV. The critical temperature in this case is Tc — 40 MeV. Left panel: mass of the condensed neutral 
kaon M2 (upper blue solid curve), and SM2 oc </> (lower blue solid curve) as well as mass of the uncondensed charged kaon Mi 
(red solid curve) in units of the chemical potential of both kaons, here chosen to be equal (and shown as the black dashed line). 
Right panel: energy gaps of the modes ef and £2 (the energy gaps e~ are larger and not shown in the figure, cf. similar plot in 
Fig-El). The condensed mode (blue curve) is gapless below the critical temperature, in accordance with the Goldstone theorem. 




FIG. 6: (Color online) Critical temperature and nonzero-temperature phase diagram for bare kaon masses as in Fig.O mj(+ — 5 
MeV, m^fo = 4 MeV. Left panel: Tc in units of MeV as a function of /i^o for fixed charged kaon chemical potentials = 
0.2, 3.2, 6.2, 9.2 MeV from top to bottom. Right panel: nonzero-temperature phase diagram in the /i^+-/ijfO-plane. From left to 
right and bottom to top (in the order black, red, blue, green, violet) the phase transition lines correspond to T = 0, 40, 50, 60, 70 
MeV. 



We now solve Eqs. ([7^ numerically from T = up to T = 2Tc. In Fig. [5] we present the self-consistent kaon masses 
and energy gaps for bare masses mi — 5 MeV, 7712 = 4 MeV, and chemical potentials /ii = 1^2 ~ 4.5 MeV. As can 
be seen from the phase diagram in Fig. [1] this corresponds to a condensed phase at T = 0. We find that the 
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critical temperature is T ~ 40 MeV, and the condensate at zero temperature ({>/ f-R — 0.6. Since our expansion in the 
condensate requires (j)/ fir ^ 1; this parameter region is close to the limit above which we do not trust the results 
quantitatively. The qualitative behavior of the quantities shown in Fig. [5] is easy to understand in view of the results 
of (^^ theory in the previous section. The uncondensed kaon, absent in the scalar (^^ theory, exhibits almost constant 
masses and energy gaps below Tc. Its energy gap is nonzero since we have chosen the charged and neutral kaon masses 
to be different, mi ^ m2. As mentioned above, if isospin is an exact symmetry, then both energy gaps f-^ik ~ 0) and 
ej(fc = 0) vanish. At T = 0, this is obvious from Eqs. ((8T|) . At nonzero T, we have confirmed numerically that in this 
symmetric case both energy gaps vanish for all temperatures below Tc. 

We finally extract the critical temperature and extend the phase diagram from Fig. [1] to nonzero temperatures. In 
Fig. [S] we show the critical temperature as a function of the neutral charged kaon chemical potential for several fixed 
charged kaon chemical potentials. Two of the curves start at the neutral kaon mass. This is clear since there is no 
condensation for chemical potentials below this mass. The other two curves for > start at a larger value, 

given by the relation (|24p . Below this critical value there is a condensate, but no condensate (the plot only 
shows the critical temperature). The jump of Tc at this critical chemical potential indicates a first order phase 
transition. 

The right panel shows the phase transition lines for both /v ° and A' + condensates. Here, the case of a charged kaon 
condensate is obtained from the above formulas simply by exchanging the indices 1 and 2. The phase transition lines 
are shown for temperatures between T = and T = 70 MeV. Due to the melting of the kaon condensates the region 
of the CFL phase becomes larger with increasing temperature. In contrast to Fig. [1] "CFL" does not mean absence 
of kaons because at nonzero temperature a thermal kaon population is present. The phase transition lines indicate 
transitions from a kaon-condensed state to a state without kaon condensation or from a charged kaon condensate to 
a neutral kaon condensate. Transitions from one condensate to another are of first order (given by the solid line). It 
is interesting that the transition from a charged kaon condensate to the CFL phase is always of second order (second 
order phase transitions are given by dashed lines). The transition from the neutral kaon condensate to the CFL phase, 
however, can be either of first or second order. It is of first order where this transition occurs on the diagonal line - 
more precisely, the line given by Eq. (|24p . Otherwise, it is of second order. This is equivalent to the statement that 
for instance the two (green) lines corresponding to T = 60 MeV do not meet the line /ii(/i2 — "^1) = ^2(^1 — rn\) at 
the same point. The do so only for equal meson masses mx+ = m^^o. Hence, the asymmetry in the order of phase 
transitions is created by the larger mass of the charged kaon. 

Since we have not imposed any neutrality constraint in this section (wc will do so in the next section), most of the 
states in the phase diagram have an overall charge. Electric neutrality fixes the charge chemical potential hq. Thus, 
if one makes use of the kaon chemical potentials in Eq. (|20p . the neutrality condition fixes, for any T, a certain point 
in the phase diagram (for a fixed quark chemical potential /i). 

Here, ignoring the neutrality constraint, we obtain an estimate for the critical temperature for arbitrary values of 
the kaon chemical potentials. The critical temperature rises rapidly with the chemical potential. In our numerical 
example {mK+ = 5 MeV, nixo = 4 MeV), when the chemical potential is only about 10% larger than the boson mass 
{fJ-KO = 4.5 MeV), the critical temperature has already risen to 40 MeV. If both kaon chemical potentials are about 
twice the kaon mass (/i/f+ ~ — 8 MeV) then the critical temperature is about 70 MeV. Of course, as mentioned 
above, the reliability of the effective theory at these large temperatures is questionable since fermionic degrees of 
freedom start to become important or the CFL phase itself has already melted. 

V. EFFECTS OF ELECTRIC NEUTRALITY 

In the previous sections we have not included any neutrality constraint. Some of the systems we have considered 
have a nonzero overall electric charge. For instance, in the zero-temperature phase diagram in Fig. [1] all states with 
if + condensate are positively charged, while all others (the if^-condensed states as well as the pure CFL state) are 
neutral. In the case of nonzero temperatures, all states we have considered are charged due to a thermal population of 
charged kaons. In this section, we impose the constraint of overall electric neutrality and include contributions from 
electrons (and positrons). Our main interest is the effect of the neutrality constraint on the critical temperature. 

A. Kaon densities 

In order to include electrical neutrality wc have to compute the kaon number densities. In fact, wc only need the 
density for the charged kaons, but in the following wc shall provide the densities of both charged and neutral kaons. 
To this end, wc have to compute the (negative of the) derivative of the effective action Vcs with respect to the kaon 
chemical potentials /xi, /X2 and evaluate the result at the stationary point. The effective potential depends on the 
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chemical potentials explicitly, as well as implicitly through cf) and S. However, at the stationary point the derivatives 
of Ves with respect to </> and S vanish by definition. Therefore, the densities can be computed from taking the explicit 
derivative of Vcff with respect to /ii, fi2- From the effective potential given in Eq. (j43|) we thus conclude 



dU 



-Tr 



(85) 



Here, U is obtained from Eq. (|62p by setting 0i = and defining = 02. The inverse tree-level propagator Sq^ and 
the full propagator S are defined in Eqs. ((64|) and ((68|) . respectively. Note that the form of the density ([85|) does not 
depend on the specific approximation of the set of 2PI diagrams V2 ■ We have only used the fact that V2 depends on 
the chemical potentials implicitly through (/) and S, but not explicitly. 
Now we can straightforwardly compute the densities. For ni we use 



iTr 
2 



K 



(86) 



where, in the second step, we have performed the Matsubara sum. For 712 we need 
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(87) 



We have made use of the excitation energies in Eq. ([70)1 with SAIi = 0. We have defined = y '^Ki^tiE^) + '^i i^)]^ 
as in Eq. ([80|) . and J,; = ■^Y^Ki^t i-^) ~ '^i^(^)]- For the explicit forms of It and Ji after performing the Matsubara 
sum see Eqs. (|B4p and (|B5|) . Now, dropping the vacuum contribution, adding the term from the potential U, and 
neglecting J2 as well as the terms proportional to 6M2 in the square roots on the right-hand side of Eq. ([87|) . we 
arrive at 



Til 



n2 



E' 



A^20^ 1 



(2^^^^2)-Y 1^^/1 + 2^/2 



(88a) 
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B. Including electric neutrality 



The electric charge density is 



/cP k 
^^/F(eci), 



where we included the electron contribution with the Fermi distribution function fpix) 
excitation energies 



(89) 

l/(e^/'^ + 1) and the 



(90) 



19 



with the electron mass rrici. Here, ejjj corresponds to positrons, while e~j corresponds to electrons; in other words, /ig 
is the chemical potential for positive electric charge, as introduced in Eq. (H^. We assume that there is no lepton 
chemical potential. Now we can solve the self-consistency equations ([79|) together with the neutrality condition 

= (91) 

for the variables (p, Mi, M2, and /xg. In contrast to the previous subsections, where we have fixed fix ~ /ig + m^/ (2/i) 
and ^2 — ■we now solely fix the parameter m^/(2/x). Then, the sclf-consistently determined charge chemical 

potential /ig determines /ii, which hence becomes a function of temperature. 




20 40 60 80 100 20 40 60 80 100 

T[MeV] T[MeV] 



FIG. 7: (Color online) Effect of the neutrality condition and the presence of electrons on the self-consistent masses and the 
critical temperature. The left panel compares the masses with neutrality condition (solid) with the ones without neutrahty 
condition (dashed-dotted). The dashed-dotted lines are taken from Fig.[S] As in Fig. [3 for both sets of curves the bare masses 
are mi = 5 MeV, 7712 — 4 MeV, and the chemical potential of the condensed neutral kaon is fi2 = 4.5 MeV. The chemical 
potential of the uncondensed charged kaon is chosen to be identical to 112 sd. T — 0. The right panel shows the negative of the 
charge chemical potential pQ, and the charged kaon chemical potential ni — /iq + /i2. 



In Fig. [7] we show the self-consistent masses as well as the chemical potentials for the neutral system. We compare 
this to the charged system that has the same kaon chemical potentials at T = 0. We see that, compared to that 
system, neutrality leads to a larger critical temperature, i.e., the kaon condensate is more robust. One can understand 
this result intuitively from the phase diagram in Fig. [S] The phase transition lines in the upper left part of this 
phase diagram imply that the K'^ condensate becomes more robust for decreasing and fixed /ii<-o, i.e., following 
a horizontal line from right to left. The reason for this is that the phase transition lines are not horizontal. In the 
neutral system, this is exactly how the chemical potentials behave as a function of temperature: fiK° is fixed while 
fiji+ decreases because of an increasing —fig- Therefore, the critical temperature is larger than that of the system 
where both fx^o and (1]^+ are fixed (at the values that the corresponding quantities in the neutral system assume at 
zero temperature). For the given parameters, we find Tc ~ 47 MeV for the neutral system, while Tc ~ 40 MeV for 
the charged system. 



VI. CONCLUSIONS AND OUTLOOK 



We have investigated kaon condensation in the CFL phase at nonzero temperatures. This study provides the basis 
for future calculations of transport properties of possible kaon-condensed CFL phases in the interior of compact stars. 
We have used the effective theory which, analogous to usual chiral perturbation theory, contains the Goldstone modes 
as dynamic degrees of freedom. 

The nonzero-temperature study of this bosonic system requires a self-consistent scheme in order to compute the 
temperature-dependent masses. Simple approximations fail since in this case the melting of the condensate would lead 
to unphysical negative excitation energies. We have thus employed the CJT formalism, which provides self-consistent 
equations through stationarity conditions for the two-particle irreducible effective action. We have treated the kaon 
chemical potentials as free, but fixed, parameters and have determined the kaon masses as well as the condensate 
self-consistently for arbitrary temperatures. This has been done in an approximation of the effective theory which 
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reduces the Lagrangian to a form not unlike usual ip'^ theory for a scalar field. The main differences to ip'^ theory are 
the coupling constants, which turn out to be combinations of the kaon chemical potentials and bare masses and the 
presence of a second, uncondensed mode. 

Our main result is the critical temperature of condensation, for which we have derived an analytic approximation. 
This high-temperature approximation reproduces the full numerical results to high accuracy. For likely values of the 
bare kaon masses of a few MeV and kaon chemical potentials slightly larger, the critical temperature is at least of 
the order of 20 MeV. Kaon chemical potentials of roughly twice the kaon masses lead to larger critical temperatures 
of approximately 80 MeV. This means that a (neutral) kaon condensate persists up to temperatures where fermionic 
excitations become important and the effective theory is no longer valid. We conclude that, if the kaon chemical 
potentials and bare masses favor kaon condensation at T = 0, the condensate will be present in a compact star which 
typically exhibits temperatures smaller than 1 MeV. 

Our results can be used for the calculation of several physical properties of the kaon-condensed CFL phase, for 
instance the specific heat. One can also compute neutrino emissivities, including both condensed and thermal kaons 
on a consistent footing. Both properties play an important role for the cooling properties of the star. So far, these 
calculations have been done without taking into account temperature effects in the Goldstone boson excitations [20| . 
Furthermore, one can use our results to compute the bulk viscosity of CFL quark matter. This has been done for the 
case where kaons do not condense [l^j, and for the superfluid mode alone [3]. Furthermore, one might include the 
effect of neutrino trapping into our calculations, which possibly leads to coexistence of two meson condensates [29| . 
Finally, it would be interesting to include higher orders in the meson fields into the calculation (at zero temperature, 
it is easy to include all orders). 
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APPENDIX A: ANSATZ FOR SELF-CONSISTENT BOSON PROPAGATOR 

In this appendix, we prove that the ansatz for the inverse propagator S^^ in Eq. (|47|) is the most general possibility. 
(The proof is presented for ip'^ theory, but the case of kaons in CFL is, within our approximations, completely 
analogous.) Let us derive the stationarity conditions without specifying any of the elements of the propagator. We 
denote the elements of the inverse propagator 5*^^ as follows, 

^ ^'^^-y S-\K) S-,\K) ) ^ - detS-HK) [ ~S^,\K) S^,\K) ) ' ^^^^ 

Then, the first stationarity equation (|46ap reads 

Q 

where we used the tree- level propagator (HH). The Dyson-Schwinger equation (j46b|) yields for the diagonal elements 

The first of these equations, Eq. (|A3a|) . can be inserted into Eq. (|A2p to obtain S^i{K) = —K^ + 2A(/)^. This can be 
re-parametrized without loss of generality as 

{K) = -K^ + Ml - I? . (A4) 

(Yielding the condition — /.j^ = 2A(^^ as one of the self-consistency equations.) We may now subtract Eqs. (|A3a[) 
and (jA3bp from each other. 



(A3a) 
(A3b) 
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Since the self-energy is momentum independent (of. its diagram in Fig. [2]) the right-hand side of this equation has 
to be momentum independent. Hence, the if-dependence on the left-hand side of Eq. (|A5[) has to cancel out. In 
other words, S^i^K) and S22{K) must have the same momentum dependence. Choosing a parametrization of the 
momentum independent part of S22 [K) we thus conclude from the form of S'f/ [K) in Eq. (|A4p 



S-2\K) = -K^ + Ml-pi\ 



(A6) 



Having fixed the form of the diagonal elements, we turn to the off-diagonal elements of the Dyson-Schwinger 
equation. After adding and subtracting these equations to and from each other, we obtain 



S^2\K) + S2,\K) 



Q 



det5'-i(g) 



(A7a) 
(A7b) 



Again, we use the fact that the right-hand side of Eq. (jA7ap and thus also the left-hand side is momentum independent. 
Consequently, the momentum-dependent parts of 8^2 i-^) ^^^d S2i{K) must have the same absolute value and opposite 
sign. On account of Eq. (jA7bp these parts must be zL2ikQiJ,. Also on account of Eq. (|A7bp . the momentum- independent 
part of both off-diagonal elements must be the same. Hence, we can write 



Inserting this back into Eq. (|A7a[) yields 



= c 



V ^ det^-i(Q) 



(A8) 



(A9) 



The expression in square brackets does in general not vanish (note that the momentum sum vanishes for T = 0). 
Hence, c = 0. In conclusion, Eqs. (jA4[) . (|A6[) . and (|A8[) with c = define the most general ansatz, as used in the 
main part of the paper. Note that this proof is specific for the case 0^0, i.e., it is valid only below the critical 
temperature. However, above the critical temperature, very similar arguments apply and one obtains the same ansatz 



with the additional constraint M^ 



M_ 



APPENDIX B: MATSUBARA SUM 



In this appendix we perform the Matsubara sum that is needed to derive Eqs. (|57p for ip^ theory. The Matsubara 
sums for kaons are analogous. Wc rewrite the sum over bosonic Matsubara frequencies go = —2imTT into a contour 
integral in the complex energy plane. 



tVs-^ = I dLoS^{Lo)coih — 

^ 2TTiJc 2T 

90 •^'^ 



c 



/ dc.5±Hcoth-^, (Bl) 
2m J_,^+r, 2T 



where S^{ijj) is the propagator from Eq. ([51 



- (6+)2][.;2„(,-)21 



The closed contour C encircles the poles of coth[uj / {2T)\ on the imaginary axis, such that no poles of S^{u)) are 
in the enclosed region. Then, this contour can be deformed to obtain the second line in Eq. (|Bip . where > is 
infinitesimally small and where the symmetry of S^{uj) and the antisymmetry of coth[a'/(2T)] in to has been used. 
Next, we can close the resulting contour in the positive half plane. This includes the positive poles of 15*^(0;) and we 
obtain with the residue theorem 
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Consequently, 



r^[^+(Q) + ^-(Q)]^l^l 1-e coth^g,, (B4) 

2 ^ //i,,2 Z72 I A/\ f4 / 2J 



and 



1 A A 1 ,e 

t5^[5+(Q) - S-{Q)] = - J2 ^ . (B5) 



With coth[a;/(2T)] = 1 + 2/(x) and dropping the vacuum contribution one arrives at Eqs. ((57)) with the definitions 

APPENDIX C: SYMMETRY NON-RESTORATION? 

The analytical approximation of the critical temperature in Eq. (|83p seems to allow for negative expressions for T^. 
In this appendix, we discuss possible implications of this apparently unphysical result. 

The approximation ([55]) was obtained under the assumption 112 > m-2 (condensation of the neutral kaons). Conse- 
quently, a + 2/32 < leads to < 0. This absence of a real critical temperature seems to suggest that the condensate 
"refuses" to melt and persists for arbitrarily large temperatures. In fact, this possibility has been discussed in the liter- 
ature for two-field models which are equivalent to our expansion (j62p and has been termed "symmetry non-restoration" 
(SNR), see Ref. [45[. A related phenomenon can already be seen from the original self-consistency equations (|79p . 
Employing the high-temperature approximation ([82)) . we see that the self-consistent mass Mi of the uncondensed 
mode decreases linearly in T if a -I- 2/?i < 0. This suggest the onset of condensation at some large temperature, 
which has been termed "inverse symmetry breaking" (ISB) [i^. Also the self-consistent mass of the condensed mode 
behaves unexpectedly for certain parameters. From Eq. (|79cp we see that the temperature-dependent term propor- 
tional to Ii ~ I2 decreases linearly in T if a -I- 2/32 < 0. This suggests that, for large temperatures, the condensate (f> 
has to increase with increasing T in order to avoid a negative mass squared. On the other hand, we know that our 
approximation breaks down for large condensates because we have employed an expansion in 0, see comments below 
Eq. (I62p . Therefore, in our specific context, the exotic possibilities of SNR and/or ISB might be an artifact of the 
approximation, and we leave it to future studies to decide whether they remain a viable option in the full effective 
theory. 

Instead, we simply identify the parameter regions where a + 2/3^ < 0. In Fig. [5] we plot the (dashed-dotted) lines 
a + 2/3i = and a + 2/32 = in the //i-//2-pha'Se diagram. These lines have been obtained with the help of the 
definition The shaded areas to the left of (and below) these curves correspond to a -I- 2/3^ < 0. 

Regions shaded in light grey indicate ISB, i.e., a decreasing self-consistent mass of an uncondensed mode, suggest- 
ing high-temperature condensation for this mode when the mass is sufficiently small. This does not occur for the 
parameters chosen in Fig. [5l whereas the phase diagram in Fig. [6] contains such parameter regions. However, the 
smallness of the coupling constants a and (3i prevents the self-consistent mass from becoming close to zero. Hence, a 
possible onset of condensation would appear at temperatures far beyond the validity of the effective theory, rendering 
the question of ISB for the given parameters purely theoretical. 

The region shaded in dark grey corresponds to a negative critical temperature squared, i.e., possibly SNR. The 
occurrence of this region requires a sufficiently large difference in the bare kaon masses ml > 12ni^. The line 
a + 2/32 = intersects the line /i2 = m2 at the point {^1,^2) = {\/ ml — Sm^ — 3m2, 771,2). Thus, this region is not 
present in the phase diagram in Fig. [SI where the kaon masses are chosen to be sufficiently close together. Because 
of large uncertainties in these masses a large difference in kaon masses is not excluded. Hence we cannot exclude the 
possibility of exotic phenomena, and future studies going beyond our simple expansion in the fields are required. 
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